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Abstract 


Employing Biot’s theory of poroelasticity, three-dimensional vibrations in a poroelastic solid that is subjected to static 
stresses is investigated. By considering second-order coupling between stress and strain, pertinent governing equations 
are derived. A frequency equation is obtained in the case of static uniaxial stress or strain. Phase velocity against static 
uniaxial stress is computed in the case of two poroelastic solids and results are presented graphically. 
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|. Introduction 


The theory of poroelastic media originates from the 
requirements of specific problems of geophysics such 
as the problems of propagation of seismic waves. 
Detailed analysis of wave propagation in elastic 
media is given in the papers by Miklowitz (1960) and 
Kolsky (1963). When the body forces act on a solid, the 
state of the solid which is at rest gives static stress and 
strains. When the large static stress acts on a solid, the 
effect of body forces can be neglected (Mott, 1971). The 
static stress caused by large, externally impressed sur- 
face forces is assumed to give effects much larger than 
the effect of gravity. The examples of externally 
impressed forces are ferroelectric, piezoelectric mater- 
ials which have large static stress produced in them by 
some external agency, such as a battery or coil (Mott, 
1971). The elastic motion of an isotropic medium in the 
presence of body forces and static stresses is investi- 
gated (Mott, 1971). In the said paper, Mott derived 
the equations of motion in an elastic medium in the 
presence of body forces and static stresses. Further, 
elastic waveguide propagation in an infinite isotropic 
solid cylinder that is subjected to a static axial stress 
and strain is given in the paper (Mott, 1973). In the said 
paper, the effect of static axial stress and strain upon 
the velocity of the lowest-order flexural mode in solid 
circular cylinders is discussed and it has been proved 
that flexural waves in cylinders and transverse waves in 
stretched strings are of the same nature. With regard to 
poroelastic solids, the governing three-dimensional 
equations of poroelasticity are developed in the 


frequency domain (Biot, 1956). Various aspects of 
wave propagation from the perspective of geophysics 
and seismology are discussed (Tolstoy, 1973) in the 
framework of Biots theory of poroelastic solids. 
Consideration can be given to flexural vibrations in 
poroelastic solids such as beams, plates, and shells 
which are three-dimensional in nature. Flexural vibra- 
tions of  poroelastic plates are investigated 
(Theodorakopoudos and Beskos, 1994). Edge waves 
in poroelastic plate under plane stress conditions are 
studied (Reddy Perati and Tajuddin, 2003). In the 
said paper, the governing equations of plane stress 
problems in poroelastic solids are formulated. 
Flexural wave propagation in coated poroelastic cylin- 
ders with reference to fretting fatigue is investigated 
(Shah, 2011). Flexural vibrations of poroelastic circular 
cylindrical shells immersed in an acoustic medium are 
studied (Shah and Tajuddin, 2010). There is much focus 
on wave propagation in poroelastic cylinders and por- 
oelastic half spaces (Tajuddin and Ahmed, 1991; Reddy 
Perati and Tajuddin, 2000; Reddy and Tajuddin, 2000, 
2010; Shah and Tajuddin, 2010; Nageswara et al., 2011; 
Rajitha et al., 2012) in the framework of Biot’s theory 
(Biot, 1956). However, the above listed publications do 


Department of Mathematics, Kakatiya University, India 


Received: 29 April 2013; accepted: 21 September 2013 


Corresponding author: 

Malla Reddy Perati, Department of Mathematics, Kakatiya University, 
Warangal-506009, India. 

Email: mperati@yahoo.com 


Downloaded from jvc.sagepub.com at Bobst Library, New York University on April 24, 2015 


Journal of Vibration and Control 0(0) 


not contain investigation of flexural vibrations of por- 
oelastic solids in the presence of static stresses. 
Therefore, in the present paper, this is investigated in 
the framework of Biot’s theory (Biot, 1956). Here the 
phrase *poroelastic solids" in the title stands for por- 
oelastic bars, poroelastic plates, and poroelastic shell 
wherein flexural vibrations occur. The pertinent gov- 
erning equations are not readily available in the litera- 
ture. Thus, the said equations are derived from the 
corresponding equations of the general three-dimen- 
sional problems. Non-dimensional phase velocity is 
computed as a function of static uniaxial stress and 
wavenumber. 

The rest of the paper is organized as follows. In 
Section 2, formulations of governing equations are 
given. Equations of motion in the presence of static 
uniaxial stress are given in Section 3. Numerical results 
are described in Section 4. Finally, conclusions are 
given in Section 5. 


2. Governing equations 


Let ufu, v, w) and U(U, V, W) be the displacements of 
solid and fluid, respectively, in the X, Y, Z directions. 
Because of the second-order coupling between the stres- 
ses and strains, it is shown in the paper (Mott, 1971) 
that the effective stresses must be inserted in place of 
the usual stresses in the constitutive relations and the 
equations of motion. Whatever may be the boundary 
conditions, this effective stress must replace the usual 
stress at the boundary. The effective stress-strain rela- 
tions are (Mott, 1971) 


; ðu ðu 
O = Oxx — Oxy 5 XZ p 
j av ðv 
Oxy = 0 Oxx ax — Oxz dz 
j ow ow 
Ox; = Oxz — Oxx ax Oxy ay 
j ov ov 
Oy, = Oyy — Oyx dx — Oyz az 
; Qu Qu 
Opx = Oyx — Oyy rm — Oy: 3z (1) 
; ow ow 
Oyz = Oyz — Oyx ax Oyy ay 
j ou ou 
Ox Ozx Ozy ay — Oz ES 
; ov ov 
97, = 05 Oma, ZZ 
a ox Oz 
; ow ow 
o = Ozz — Og m 
K ox oy 


As there are no shear components for the fluid pres- 
sure, there will not be any effective fluid pressure as that 
of solid. The expression for the fluid pressure remains 
the same. In all the above, usual stresses o; and fluid 
pressure s are (Biot, 1956) 


Oi; = 2Nej + (Ae + Q&)ó; 
s — Qe 4 Re. 


(i,j = Xy, z), (2) 


In equation (2), ej?s are strain displacements, ô; is 
the well-known Kronecker delta function, e and € are 
the dilatations of solid and fluid, respectively, and 
A,N,Q,R are all poroelastic constants. Substitution 
of effective stresses for usual stresses, the equations of 
motion appear in the following manner: 


; pac da gd j [424 ., 9w 
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de de 
l e — (au + px), 
de de oP 
95 ay gg P" + 9n P), 
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(3) 


In equation (3), o;(i = j= x,y,z) are the effective 
stresses, pj are mass coefficients such that the sums 
Put i» and p2 + pn represents the mass of solid 
and fluid per unit volume of bulk material following 
(Biot, 1956). The coefficient pi» is a mass coupling par- 
ameter between solid and fluid phases. Furthermore, 
the mass parameters obey the inequalities oj, > 0, 
p2 7 0, p22 > Ü (911922 — p72) 2 0. and Fy, F, and 
F. are the components of the body force vector F. 
Equations (3) together with equations (1) are to be 
satisfied at every interior point of the body, and on 
the surface of the body. Appropriate boundary condi- 
tions are to be imposed so that problem can be solved 
completely. Thus, in the presence of static stress, the 
solution of a dynamical problem in poroelastic solids 
can be completely determined. By neglecting the liquid 
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effects, that is, setting p12 =0, pr2r=90, p1 =P, 
A-— g- =i, N= p, the results of purely elastic solid 
(Mott, 1971) are recovered as a special case. In the 
above, 4, u, and p are Lame constants, and density 
of purely elastic solid, respectively. To construct the 
relationships between the time variant and the static 
poroelastic quantities, the particle displacements 
and the stresses are written in the following form 
(Mott, 1971): 

u = XO n(x, y, Z, Ont) = uo d- uid uc, 

v=} noy z ox) = Yo vi vas 

w= Y Walx, y, Z, ext) = Wo + wi Wa ss, 

U 2 UG yz ost) = Uo + Ul + Us. (4) 

V= D V(X, y, Z, Ont) = Vo + Vi + V» + er) 

W — Y W(x, y,2, Ont) = Wo Wi t+ Wat, 

Ojj = SCC Ey yt) = Oij + Oij + Oij, Tee 


where o, is the nth angular frequency. Using equations 
(1), (2) and (4) in the equations (3), there follows the 
static equations (n = 0) given below: 
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Similarly, the harmonic equations (n= 1) are 
obtained, which are: 
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E HL and Fy, Fy Fq are the body forces which 
have the frequency w. Equations (5) and (6) are static 
equations, and first harmonic equations, respectively. 


3. Equations of motion in presence of 
static uniaxial stress 


Assume that the static uniaxial applied stress is acting 
in the direction of z axis. Then there follows 


ijo = E ij; = 0, (i z Jj) 
Oxxo = Oyy = 0, m 
Fo = Fy, = Fa =0, 


Oz- is the applied uniform static uniaxial stress. These 
conditions are applicable for a bar with a uniform cross 
section and with its side surface stress free. Substituting 
equations (7) in equations (5), it can be seen that equa- 
tions (5) are automatically satisfied and the static 
strains are found to be ex, = eyy = —Vezz. 
Substituting these equations in the first equation of 
equation (2), the following equation is obtained: 


e = mw — Oe (8) 


In the above, Ü — d a; is Poisson ratio, and 
NCP4-24) (AN) 
Y =z is Young's modulus. 


Substituting equation (7) and equation (8) in equa- 
tion (6), the equations of motion for the first harmonic 
terms are obtained and are: 
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In equations (9), the suffix one on displacement com- 
ponents and the stresses has been dropped, and 
F = F,, = F., = 0. Due to the assumption of infinite- 
simal deformation in the linear theory of elasticity, the 
product terms can be neglected and the equations of 
motion in this case are obtained as follows: 


2 
NV?u + (4 - NÈ A ETO y 
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In all the above, C1, C», C, C4, Cs, Cg are arbitrary (Qki — 120) Ci + QU Ca + kis C3) 
constants, j is the complex unity, and k;(i = 1,2,3) is + (RKj — p20") Ca + Rl kiko Cs + kiks Co) = 0, 
the wavenumber in the ; th direction such that the QkikaCi + (ok _ pi”) C; + Qkika Cs 
wavenumber k = J/ Kj + k3 + kj. Substituting equation — | Rkik>C4 + (RK — pae?) Cs + Rkaks Cg = 0, 
(11) in equation (10), the equations of motion in terms OlkiksC + kl C3) + (Qi _ p120?) C 
of displacements are as follows: nes 
+ R(kik3C4 + kok3 C3) + (Rk; — pnw )€s = 0. (12) 
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4. Numerical results 


For the numerical work, the propagation is 
considered along z direction. In this case, kų = 
ky =0, and equations (12) reduce to the following 
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Table |. Material parameters. 


Material 
parameters al a a3 a4 dı d; d; Z 
Material- I 0.843 0.65 0.28 0.234 0.901 —0.001 0.101 3.851 
Material-2 0.96 0.006 0.028 0.412 0.887 0 0.123 2.129 
matrix form: 
1] 
Au 0 0 Ci > 6 
0 An 0 Cy | = 0, (13) $5 
0 0 43} LCs $4 
33 
where : 2 
21 
PA = 
A Nk3 k30:2, ) — NV,?o”, 5 
11 ( 3 (P +24) 9 4 s g ? 
z Oo 93 05 07 09 13 13 £$ 
PA Static uniaxial stress 
2 2 =} 
A22 = (wa (P n 2A) Hos. NV, e, 
Figure l. Variation of non-dimensional phase velocity with the 
A3 = | Pk + EA ka; pio static uniaxial stress 
33 3 N(P + 2A) 3° zzo 11 » 
k3( 1 a 2 i ; . , 
— | Ok; NP + 4 45 p120 In equation (15), c is phase velocity, and m is the 
> : non-dimensional phase velocity. Employing these 
T Oks — pw non-dimensional quantities in the frequency equation, 
Rk? — 02202 ] an implicit relation between non-dimensional phase vel- 
ocity (m), static uniaxial stress and wavenumber is 
and Vz, = —^PJ?? — is the shear wave velocity. obtained. For the numerical process, two types of 


P11 9227 £15 à . 
For a nontrivial solution, the determinant of coeffi- 
cients matrix is zero. Accordingly, the following fre- 
quency equation can be obtained: 


Au 0 0 
0 A4» 0/20. (14) 
0 0 A3 


The frequency equation is investigated by introdu- 
cing the non-dimensional quantities given below: 
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vila 


e 2 
p OF 


(15) 


materials are considered, namely, material-1 is made 
up of sandstone saturated with kerosene (Yew and 
Jogi, 1976), while material-2 is sandstone saturated 
with water (Fatt,1959). The physical parameters of 
the said materials pertaining to equation (15) are 
given in Table 1. Static uniaxial stresses under consid- 
eration are external surface forces which can vary, 
hence phase velocity values are computed against 
static uniaxial stress for a fixed wavenumber. The 
values are computed using the bisection method imple- 
mented in MATLAB and the results are depicted in 
Figure 1. Figure 1 shows the plots of non-dimensional 
phase velocity against the static uniaxial stress for fixed 
wavenumber in the case of material-1 and material-2. 
From this figure it is clear that material-1 value is 
greater than that of material-2 for the same value of 
applied static stress. Both the materials are sandstone 
related and differ in only the fluid part. Hence, it can be 
inferred that the fluid part is causing above discrep- 
ancy. When the static stress increases, non-dimensional 
phase velocity values increase for both materials. 
Moreover, for fixed static uniaxial stress, the non- 
dimensional phase velocity value is computed against 
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wavenumber, and the values found to be the same for 
all values of wavenumber. For example, when the static 
stress is 0.3, the non-dimensional phase velocity is 
4.8047 for all the values of the wavenumber in the 
range 1—20. Therefore, in this case, it may be concluded 
that phase velocity is independent of wavenumber. 


5. Conclusions 


Employing Biot’s theory, flexural vibrations of poroe- 
lastic solids in the presence of static stress are investi- 
gated (Biot, 1956). Pertinent constitutive relations and 
the equations of motion are derived. Non-dimensional 
phase velocity against static uniaxial stress for fixed 
wavenumber is computed for two types of poroelastic 
solids. Static uniaxial stresses under consideration are 
external surface forces which can vary. Therefore, the 
phase velocity is computed against static uniaxial stress. 
From the results, it is clear that material-1 value is 
much greater than that of material-2 for the same 
value of applied static stress. Both the materials are 
sandstone related and differ in only the fluid part. 
Hence, it can be inferred that the fluid part is causing 
the above discrepancy. When the uniaxial static stress is 
fixed, it is concluded that the phase velocity is inde- 
pendent of wavenumber. Similar analysis can be made 
for different poroelastic solids if the values of their por- 
oelastic constants are available. 
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